LA-UR- 12-23703 



A Dispersive Treatment of K l4 Decays 

Gilberto Colangelo 1 , Emilie Passemar 2 , and Peter Staffer 1 a 

1 Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics, University of 
Bern, Sidlerstrasse 5, CH-3012 Bern, Switzerland 

2 Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA 



Abstract K H are for several reasons an especially interesting decay channel of K mesons: 
K e4 decays allow an accurate measurement of a combination of 5 -wave nn scattering 
lengths, one form factor of the decay is connected to the chiral anomaly and the decay is 
the best source for the determination of some low energy constants of ChPT. We present 
a dispersive approach to Km decays, which takes rescattering effects fully into account. 
Some fits to NA48/2 and E865 measurements and results of the matching to ChPT are 
shown. 

1 Motivation 

Kn, the semileptonic decay of a kaon into two pions and a lepton-neutrino pair, plays a crucial role in 
the context of low energy hadron physics, because it provides almost unique information about some 
of the 0(p 4 ) low energy constants (LECs) of Chiral Perturbation Theory (ChPT), the effective low 
energy theory of QCD. The physical region of K(\ starts already at the mi threshold, thus it happens at 
lower energies than e.g. elastic Kn scattering. Since ChPT is an expansion in the masses and momenta, 
it is expected to converge better at lower energies. Therefore Km is a particularly interesting process 
to study. 

Besides, as the hadronic final state contains two pions, K^ is also one of the best sources of 
information on the nn scattering lengths a° Q and (T). 

On the experimental side, we are confronted with impressive precision from high statistics meas- 
urements. During the last decade, the process has been measured in the E865 experiment at BNL O 
and in the NA48/2 experiment at CERN Q]. Very recently, the NA48/2 collaboration has published 
the results on the branching ratio and form factors of Km, based on more than a million events Q. 

Here, we present preliminary results of a new dispersive treatment of K(4 decays. Dispersion rela- 
tions are an interesting tool to treat low energy hadronic processes. They are based on the very general 
principles of analyticity and unitarity. The derivation of the dispersion relation applies a chiral power 
counting and is valid up to and including 0(p 6 ). The dispersion relation is parametrised by five subtrac- 
tion constants. As soon as these constants have been fixed, the energy dependence is fully determined 
by the dispersion relation. The presented method implements a summation of final state rescattering 
effects, thus we expect it to incorporate the most important contributions beyond 0(p 6 ). 
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2 Dispersion Relation for K, 4 Decays 
2.1 Matrix Element and Form Factors 

We consider the charged decay mode of Kn'. 

K + (p) -> n + {p x )n-{p 2 )t{p e )v t {p v ), (1) 

where I e [e,p] is either an electron or a muon. 

After integrating out the W boson, we end up with a Fermi type current-current interaction and 
the matrix element splits up into a leptonic times a hadronic part. The leptonic matrix element can 
be treated in a standard way. The hadronic matrix element exhibits the usual V - A structure of weak 
interaction. Its Lorentz structure allows us to write the two contributions as 

<;r + (/> 1 >r(j>2)|^0)|.r( / >)> = --^e^IlpfQ?, (2) 

K 

(n + (p l )n-(p 2 )\A IJ (0)\K + (p)) = -ijL (p^F + Q^G + L^R) . (3) 

In the electron mode (up to now the only one where experimental data is available), mainly one specific 
linear combination of the form factors F and G is accessible: 

PL 

Fi{s,t,u) = XF(s,t,u) + {u- t)—G(s,t,u), (4) 

where s, t and u are the usual Mandelstam variables, X = ^A l ^ 2 {M 2 ( , s, S(), PL = \(M 2 K — s - S() and 
si = (pc + Pv) 2 - A(a, b, c) — a 2 + b 2 + c 1 - 2(ab + be + ca) is the Kallen triangle function. 



2.2 Decomposition of the Form Factor 

The form factor F\ has the following analytic properties: 

- There is a right-hand branch cut in the complex s-plane, starting at the 7T7r-threshold. 

- Analogously, in the t- and M-channel, right-hand cuts start at the /T7r-threshold. 

Due to crossing, the right-hand cuts in the t- and M-channel show up in the ^-channel for negative 
values of s. The situation is analogous for the other channels. 

Based on fixed-? and fixed-M dispersion relations for the form factor, we can derive its decompos- 
ition into functions of a single variable. Such a decomposition has first been worked out for the nn 
scattering amplitude [4] and later for Kn scattering (5). 

The basic idea is to define functions that only contain the right-hand cut of each partial wave and to 
split up in this way all the discontinuities of the form factor. For instance, in the case of the ^-channel 
S -wave /o, this function looks as follows: 

M (s):=P(s) + S - f 2 ( , lmMs :\ >4 ds', (5) 
71 JamI - s- ie)s'* 

where P(s) is a subtraction polynomial. After defining similar functions that take care of the right- 
hand cuts of f\ and the S - and P- waves in the crossed channels, all the discontinuities are divided into 
functions of a single variable. 

This procedure, also known as 'reconstruction theorem', neglects on the one hand the imaginary 
parts of D- and higher waves, on the other hand the high energy tails of the dispersion integrals from A 2 
to oo. Both effects are of 0{p & ) in the chiral counting. In order to simplify substantially the dispersion 
relation, we also neglect at the present preliminary stage the dependence on Sf , which is experimentally 
small. (Otherwise, we would need to consider a coupled system for the form factors F and G.) 

Respecting isospin properties, we obtain the following decomposition of the form factor: 



2 I 2 

FiOs, t, u) = M (s) + -No(t) + -R (t) + R (u) + (u- t)M x (s) - -[t(u - s) - A Kn Aen\Ni{t), (6) 

where A K „ = M 2 K - M 2 and A Cn - s e - M 2 . 
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2.3 Integral Equations 

One should not try to solve directly the dispersion relation |5]l since it may not uniquely determine the 
solution of the problem [6]. Noting that each of the functions Mo, . . . satisfies an Omnes equation, we 
apply the solution to the inhomogeneous Omnes problem: 



Similar relations hold for the other functions. Due to Watson's final state theorem, we can identify the 
phases 5\ (I - angular momentum, / - isospin) with the phase shifts of the elastic scattering. As an 
input to our equations, we therefore need the following phase shifts: 

- <5q, 6\: elastic mi scattering [7| 

- dl' 2 , d\' 2 , 6 3 J 2 : elastic Kn scattering 15191 

The inhomogeneities in the Omnes problem are given by the differences of the functions Mo, . . . 
and the corresponding partial wave, e.g. M)(s) = fo(s) - Mq(s). These 'hat functions' contain the left- 
hand cut of the partial wave and we compute them by projecting out the partial wave of the decomposed 
form factor (|6J. E.g. Mq{s) is then given as angular averages of No, Ni, etc. 

We can now solve the dispersion relation for the form factor. We have parametrised the problem 
by the constants appearing in the subtraction polynomials (in total, only five independent subtraction 
constants are needed due to an ambiguity in the decomposition) and we use the elastic scattering phase 
shifts as inputs. The energy dependence is then fully determined by the dispersion relation. 

We face a set of coupled integral equations: The functions Mq(s), Mi(s), . . . are given as dispers- 
ive integrals involving the hat functions Mq(s), Mi(s), . . ., whereas the hat functions are themselves 
defined as angular integrals over Mo(s), M\ (s), etc. This system can be solved by iteration. The prob- 
lem is linear in the subtraction constants that have to be determined by a fit to data. 

3 Preliminary Results 

3.1 Fit to Data 

We perform a fit of the dispersion relation to both, the E865 1 2 1 and NA48/2 data sets [ 3 1 . The S -wave 
dominantly determines the three subtraction constants in Mo, the f-wave the two in M\ . Figure [T] 
shows the result of the combined fit for the two partial waves. The^ 2 /dof of this fit is 1.73. This rather 
large value is a consequence of the small uncertainties of the relative form factors in the NA48/2 data. 

3.2 Matching to ChPT 

After having fixed the subtraction constants by fitting the data, we can match our dispersive represent- 
ation to the ChPT result. In this matching procedure, we determine the LECs L[, U, and U y 

The dispersion relation allows us to choose a convenient matching point that may lie outside the 
physical region. We choose the point s = t — u = below the threshold, where ChPT should even 
converge better. In table [T] we show preliminary values of the low energy constants, resulting from 
a matching to the 0(p 4 ) ChPT result. For comparison, we also quote the values of the global fit by 
Bijnens and Jemos 1 10 1 that takes as input the Km data from NA48/2 |[T). 

We find that our fit to NA48/2 agrees with this global fit. A large part of our uncertainties stems 
from the nn scattering phase shifts and contains at present a conservative estimate of the systematics. 

Our preliminary values are subject to change as we are going to perform a matching to 0(p b ) ChPT 
and include the neglected st dependence as well as isospin breaking corrections ifTTl . 




with a new subtraction polynomial P(s) and the Omnes function 




(8) 
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6.3 




s [Mj] s [Mj] 

Figure 1. Fit of the 5- and P-wave to both E865 [2| and NA48/2 [3| data. The error bars indicate separately the 
uncorrelated part of the uncertainties and the total errors (containing also the uncertainty from the normalisation). 



Table 1. Preliminary results for the LECs (fi = 770 MeV). 





lO 3 ^ 


10 3 L!J 


10 3 L^ 


Dispersive treatment, fit to E865 J2) 


0.42 ± 0.41 


0.41 


+ 0.34 


-2.19 ± 


1.41 


Dispersive treatment, fit to NA48/2 |3] 


0.60 ± 0.29 


0.63 


±0.28 


-3.16 ± 


1.19 


Bijnens, Jemos, 'fit AH' 1101 


0.88 ± 0.09 


0.61 


±0.20 


-3.04 ± 


0.43 



3.3 Summary 

We have presented a dispersive representation of K(a decays that provides a model independent para- 
metrisation valid up to and including 0(p b ). It includes a full summation of final state rescattering 
effects. It is parametrised by subtraction constants that we fix by fitting experimental data. The match- 
ing to ChPT can be performed below the physical threshold, where ChPT should converge better. 
Hence, we expect to find more reliable values for the LECs than with a pure ChPT treatment. 
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